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Lieber Derive Anwender! 

Das Echo auf D-N-L# 1 war gewichtig, nicht so sehr 
in quantitativer als vielmehr in qualitativer Hinsicht. 
Bis jetzt (7.April) hat die DUG 128 Mitglieder, dar- 
unter viele Hochschulprofessoren, Hochschulinstitute 
und Anwender aus Industrie und Forschung. Ich 
ersuche alle Derive User - besonders jene aus Lehre, 
Forschung und/oder Industrie mit uns iiber ihre spe- 
ziellen Derive- Anwendungen in Kontakt zu treten. So 
konnen wir die Bandbreite der Moglichkeiten und 
Grenzen von Derive fur alle allmahlich erweitem. 

In dieser Ausgabe werden zwei Beitrage been- 
det, ein neuer - iiber DGL begonnen. Im weiteren 
sollen in loser Folge die Utility-Files besprochen 
werden. Vielleicht konnten Mitglieder der User 
Group spezielle Files besprechen, bzw. verbessem 
oder erganzen? 

Einen Schwerpunkt will ich noch besonders betonen: 
Mathematik ist schon, das wissen die Mathematiker 
nicht erst seit dem Blick in die faszinierende Welt der 
Fraktale und des Chaos. Auch die Schonheit eines 
Beweises, einer Formulierung, eines Algorithmus 
kann begeistem. Derive laBt dank seiner kombi- 
nierten rechnerisch-graphischen Fahigkeiten die 
Schonheit der Mathematik erkennen. So sind in diese 
Ausgabe Graphiken eingestreut, die leicht mit dem 
VECTOR-Befehl als Kurvenscharen erzeugt werden 
konnen. Ich mochte Sie auffordern, selbst in dieser 
Richtung zu experimentieren. Lassen Sie sich iiberra- 
schen und schicken Sie mir bitte Ihre "schonsten" 
Bilder - was immer das heiBen mag. Ich glaube, daB 
auch Schuler und Studenten (und Lehrer!) an diesem 
Experimentieren mit Winkel-, Exponential- und an- 
deren Funktionen, mit und ohne Parameterdarstellung 
groBen Nutzen ziehen konnen. Spielerisch lassen sich 
Eigenschaften und GesetzmaBigkeiten dieser Funk- 
tionen erkennen und spater nutzen. 

Nun noch eine Entschuldigung: aus einem tech- 
nischen Versehen wurde bei D-N-L# 1 das Erschei- 
nungs-datum mit Janner 1990 angegeben; naturlich 
war der Janner 1991 gemeint. Dies ist nun die echte 
Nummer 2 aus April 1991. 

Von einigen wenigen Mitgliedem ist bis jetzt 
nur die Anmeldung eingelangt, die wir mit Freude 
registriert haben. Leider wurde noch nicht der Mit- 
gliedsbeitrag beglichen. Wir bitten Sie, den Mit- 
gliedsbeitrag (oS 300.- fur Osterreich, bzw. DM 45.- 
als Scheck vom Ausland - sonst sind die Bankspesen 
zu hoch) zu entrichten. Der Mitgliedsbeitrag dient zur 
Ganze der Begleichung der Gestehungskosten des 
D-N-L. 

Mit den besten GriiBen 



Dear Derive User! 

D-N-L #1 has met with a fairly lively response, not 
so much in quantitative than rather in qualitative 
respect. Up to now (7 April) 128 members have 
joined the DUG. Among them there are a lot of uni- 
versity professors and university institutes, teachers 
as well as users coming from research and industry. I 
would like to ask all Derive Users, especially those 
coming from teaching, re-search and/or industry, to 
contact us with regard to their experiences made with 
their special applications of Derive. In this way we 
will be able to expand the range of Derive's possibili- 
ties and limits by and by. 

In this issue two contributions are finished and a 
new one - concerning ODEs - is started. In the fol- 
lowing issues other utility files will be reviewed. I 
wonder if some members of the DUG might discuss, 
complete or improve a special file? 

Let me now put special emphasis on the under- 
lying theme to one of the central topics of this issue: 
Mathematics is beautiful. As a matter of fact, mathe- 
maticians have come to recognize this not only since 
they started throwing glances at the fascinating world 
of fractals and chaos. There is indeed, believe it or 
not, a beauty to be found in a proof, a particular for- 
mulation or an algorithm which one may well get 
enthusiastic about. Due to Derive it is easy to per- 
ceive this beauty of maths thanks to its combined 
calculating-graphic capabilities. This is to be illus- 
trated by means of plots, which are included in this as 
well as in later issues and can easily be generated as 
families of curves by applying the VECTOR- 
command. I want to invite you to experiment yourself 
in this direction and just to be surprised by the re- 
sults. I would be glad to receive your "most beau- 
tiful" pictures, whatever this may mean. As far as I 
am concerned, I would feel that pupils, students (and 
maybe their teachers) may derive great profit from 
this working with trigonometric, exponential and 
other functions using or not the parametric form. The 
properties and regularities of these functions may be 
recognized in a playful manner and exploited for later 
use. 

Last but not least an apology: because of a tech- 
nical mistake the publishing date of D-N-L #1 was 
indicated as January 1990. What was meant, of 
course, was January 1991. This, however, is the true 
number 2 of April 1991. 

At last a reminder: a few members have sent 
their applications, which have been accepted with 
pleasure, without however settling the membership 
dues of GM 45.- (or AS 300.- in Austria). Please be 
so kind to pay your dues (send a cheque for GM 45.- 
from outside Austria because of the considerable 
banking charges) needed to cover solely D-N-L- 
production and distribution costs. Let me thank you 
in advance for remitting promptly. 

With best regards 
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The Derive-News-Letter is the 
Bulletin of the Derive-User- 
Group. It is published at least 
three times a year with a con- 
tent of 30 pages minimum. The 
goals of the D-N-L is to enable 
the exchange of experiencies 
made with Derive as well as to 
create a group to discuss the 
possibilities of new methodic 
an didactic manners in teaching 
Mathematics . 

Subscription of the D-N-L is 
restricted to Members of the 
Derive-User-Group . Membership- 
form is enclosed with the first 
issue you get. 

Editor : 

Mag. Josef Bohm 
A-3042 Wurmla 
D ’ Lust 1 
Austria 


Contributions : 

Please send all contributions 
to the above adress. Non-Eng- 
lish speakers are encouraged to 
write their contributions in 
English to confirm the interna- 
tional touch of the D-N-L . At 
the other hand non-English ar- 
ticles will be warmly welcomed, 
too. Your contributions will be 
edited but not refereed. By 
submitting articles the author 
gives his consent for reprin- 
ting in the D-N-L. The more 
contributions you will send to 
the Editor the more livelier 
and richer in content the De- 
rive-News-Letter will be. 


Preview (Contributions for the next issues): 

With ITERATES to Chaos 
How to write one's own DEMO-file 
Solving ODE's using DERIVE (Part 2) 
More graphics 


(will be published September 1991) 


Impressum: 

Medieninhaber : DERIVE User Group, A-3042 Wurmla, D'Lust 1, AUSTRIA; 
Richtung : Fachzeitschrif t ; 

Herausgeber: Mag. Josef Bohm, Adresse wie oben; 

Hers tel lung : eigene Vervielf altigung 
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Derive - User - Forum 


E . -M . Baumgardt , Reinbek 

...Ich habe Nachhilfeunterricht in alien mathematisch - naturwissenschaf t- 
lichen Fachern erteilt, solange ich aus familiaren Griinden nicht aufter Haus 
arbeiten konnte . DERIVE war mir und meinen Schiilerlnnen dabei sehr niitz- 
lich. Mit diesem Programm lieft sich auch vollig verunsicherten Jugendlichen 
zeigen, daft 

1: der Mathematik mit gesundem Menschenverstand und mittelmaftigem Ar- 

beitsaufwand beizukommen ist 

2: auch die ekelhaf testen Mathe/Physiklehrer sehr zahm werden, sobald man 

sie fachlich an ihre Grenzen drangt 

3: ein Computer zwar ein nutzliches Hilfsmittel ist, einem aber kein Stuck 

Denkarbeit abnimmt (siehe Nullstellen der Polynome vom Grad 4 und ho- 
her) . 

Im Ubrigen hat der liebe Gott die Mathematik erfunden, um die Menschheit zu 
ergotzen und nicht, um arme Schuler zu argern. 

Mein dringlichster Wunsch an DERIVE: eine Moglichkeit, Grafiken in einem 
gangigen Format abzuspeichern (PCX, TIFF, . ..) . 

Fur Herrn K. Schmidt einen besseren Wirkungsgrad ! 

D-N-L : You will find some hints for storing DERIVE-Graphics in this issue. 


StD . R . Schorn , Kaufbeuren 

Nachmittags gelesen und erschrocken, weil die Seitenzahlung chaotisch ist 
(ein Wunder der Firma X... und deren Gebrauchsanweisung) und milde emport 
uber den Brief des Herrn Dipl.-Ing. .Schmidt aus Koln. 

Die Aufterungen des Herrn Schmidt entbehren weitgehend der Sachkenntnis . Die 
Verwendung von Taschenrechnern und Computern dient uberwiegend nicht dem 
Mathematikunterricht , sondern wird bei der Berechnung von Naherungswerten 
gebraucht . Ohne Verwendung von Taschenrechnern diirfte die Beschaf tigung mit 
Trigonometrie unmoglich sein. Oder sollen wir wieder mit Logarithmentaf ein 
(habe ich vor fiber 40 Jahren machen miissen und die Freude an der Mathematik 
trotzdem nicht verloren) Oder Rechenschiebern arbeiten? Die von Herrn 
Schmidt so genannten Dummen (angeblich n-2 Schuler einer Klasse) werden al- 
so schlagartig nach dem Abitur gescheit. 0 heilige Einfalt! Da konnte man 
noch einen ganzen Roman anhangen . 


Ein Beispiel: Bei dem bekannten Beweis nach Euklid, daft es unendlich viele 
Primzahlen gibt, braucht man das Produkt der ersten n Primzahlen und ad- 
diert 1 . 


ttl : 
#2 : 
#3 : 

#4: 

#5 : 
#6 : 


PRIMS(n) := I TERRIES (NEXT_PRIME(u)„ u, 1, n) 
PRIMS (3) 

[1, 2, 3, 5] 


EUKLID(n) 


n + 1 

n ELEMENT { PR I MS (n +1), i) 

L i=l 


1 


#7: 


EUKLID{2) 

7 

EUKLID(IO) 


{und das ergibt nach FACTOR: } 







331 -571 -34231 
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#9: 


ind nach FACTOR: 


#11 : EUKLID(20) 

#12 : 557940830126698960967415391 


Dr. H.-J. Kayser, Diisseldorf 

. . . Zu den sicher wohlgemeinten Ausfuhrungen des Herrn Schmidt aus Koln, die 
"ins Grundsatzliche gingen" mochte ich mich nur ganz kurz auhern, (sonst - 
so furchte ich, wird eine unendliche Geschichte daraus) : 

Wir haben an unserer Schule (Gymnasium) sehr positive Erfahrungen mit 
DERIVE gemacht . Der Mathematik-Unterricht konnte durch den Einsatz von 
DERIVE (gelegentlich aber auch mit anderen Computer-Programmen) in Verbin- 
dung mit einem LCD-Display lebendiger und effektiver gestaltet werden. uber 
solche Erfahrungen wird demnachst in einer Schrift der Landesstelle fur 
Schule und Weiterbildung, Soest sowie in Mathematik-Zeitschrif ten, z.B. MBU 
(Mathematik Betrifft Uns, Bergmoser und Holler Verlag, 5100 Aachen) , be- 
richtet werden. Auch Lehramtsanwarter beziehen die Arbeit mit dem Computer 
und insbesondere mit dem Programm DERIVE in ihre schrif tlichen Examensar- 
beiten ein. 


Pam Bishop, CTI , University of Birmingham 

I would be delighted to accept your invitation to become founding vice 
president of the Pan-European DERIVE User Group. We have already published 
several articles about DERIVE in our own newsletters, and have set up a 
working group to share experiences of using DERIVE in higher education. 

(Pam Bishop is Information Officer of the CTI - Computers in Teaching Ini- 
tiative, University of Birmingham and University of Glasgow) 


R. Stablein, Universitat Claus thal 


Das Sof twareprodukt wird von uns regelmahig zur Uberprufung unserer Berech- 
nungen eingesetzt. Ein immer wieder auftretendes Problem stellt dabei der 
Speicherplat z des Programmes dar. So ist das Berechnen einer inversen Ma- 
trix (3x3 komplex) , das Berechnen von Eigenwerten usw. mit dem vom Programm 
zur Verfugung stehenden Speicher (DOS spezifisch nur 640 kB) auch mit lei- 
stungsf ahigeren Maschinen mit sehr viel mehr Arbeitsspeicher unseres Wis- 
sens nach unter DOS nicht moglich und somit das Ausweichen auf andere Pro- 
gramme notig (REDUCE) . 

D-N-L : See the answer from Soft Warehouse, Hawai 

The extra 286 address space and expanded memory are not very suitable for 
list oriented programs such as computer algebra - the registers are still 
16 bits and list structure requires arbitrarily tangled pointers throughout 
the entire total memory address space. Also, it is important to realize 
that going to 32-bit pointers will approximately double the memory to store 
a given expression and the size of DERIVE.EXE, making a 1.280 mb 32-bit 
system equivalent to a faster 640 kb 16Abit system. Also computer-algebra 
data space requirements often grow exponentially with problem size so that 
even a ten-fold increase of memory yields only a meagre increase in the si- 
ze of solvable problems. 

Derive will not go to virtual memory because it would then become so large 
that it would require everyone to have a 286 or better and a hard disk. 
However, you can always explicitly temporarily save expressions to files 
then remove them from the algebra window to temporarily free up space. 
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M. Lauer, Pirmasens 

Herr Lauer hat mehrere Probleme mit Derive 2.0: 

1. Was bewirkt bei ITERATE der 4. Parameter n = -1?? (Vgl . INVERSE in 
MISC.MTH) 

2. Die Programmierung von MOD in MISC.MTH sieht ja merkwurdig aus und lie- 
fert auch fehlerhafte Resultate: 

MOD (1234560, 9876) = 61.7272 ! (richtig: 60 !) 

3. Die rekursive Version von GCD (deutsches Handbuch S.225, engl . Manual 
p 214) ist nicht sehr brauchbar. 

4. Ich vermisse eine Iteration mit offenem Ende (Abbruchbedingung) : 

ITERATE (u, x, xO, b) , wobei b eine Relation ist. 

Dies laftt sich zwar mit ITERATE (IF (NOT b, u, x) , x, xO) simulieren, 
ist aber unschon. 

5. Warum stellt DERIVE nicht alle Funktionen zur Verfugung, die es selbst 
benutzt? Beispielsweise die des groLten gemeinsamen Teilers GCD, die es 
beim Kurzen von Bruchen ja braucht? 

6. Ich versuche, mit DERIVE einfache diophantische Gleichungen zu losen: 

a x + b y = c (a,b,c,x,y e Z) 

D-N-L ; Viele Fragen, daher mehrere Antworten. Hr. Lauer hat noch einige 
Bemerkungen, die nach Ruckfrage in einer der nachsten Ausgaben behandelt 
werden konnen . Nun der Reihe nach: 


Zum 

Parameter 

n = 

-1 

in 

der Funktion INVERSE aus MISC.MTH 

1: 

INVERSE 

(u. 

X) 

: = 

= ITERATE (u, x, x, -1) 

2 : 

INVERSE 

3 

(x 

+ 

1, 

x) 

3: 

(x - 1) 

1/3 




4 : 

ITERATE 

3 

(x 

+ 

1, 

1/3 

x, (x - 1) , 1) 

5: 

X 






Ich suche die Umkehrfunktion zu y = x 3 + 1 . Das ist aber jene Funktion 
g(x), die nach Anwendung der vorliegenden Iterationsvorschrif t u ebenso den 
Term x erzeugt, wie aus x der Term f (x) erzeugt wird. Das ist also gleich- 
sam der vorletzte Schritt, Oder auch der -1. Schritt. 

2. Die Funktion MOD funktioniert , wenn Sie im Exact Mode arbeiten. 

15: MOD (1234560, 9876) 

16: 60 

3. GCD (123456, 9876) funktioniert bei mir auch nicht. Interessanterweise 
geht aber GCD ( 1234560 , 987 6) schon (= 12). Ich werde dieses Phanomen an 
Soft Warehouse weitergeben (gemeinsam mit Frage 5) . 

Wenn Sie im Exact Mode arbeiten, sieht Ihr IQR(a,b) etwas bequemer aus 
(siehe Ausdruck) . Hr. Lauer benotigt IQR(a,b) als Integer Quotient and 
Remainder der Division a/b. 





p 
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17: GCD (m, n) := IF (n = 0, |m|, GCD (n, |m oo n|)) 

18: GCD (1234560, 9876) 

19: 12 

20: GCD (123456, 9876) 

21: "memory full!" 




a ) 


22: IQR (a , b) := 

FLOOR 

\ — 1 

X) 

, MOD (a, b) 


23: IQR (123456, 9876) 
24: [12, 4944] 

25: IQR (1234560, 9876) 
26: [125, 60] 


4. Ich habe Herrn Lauer gebeten, Anwendungsmoglichkeiten dieser Iteration 
beizustellen : 

Hier verdopple ich den Anfangswert 3 solange bis 100 (bzw. 1000) iiber- 
schritten werden k"nnte : 


27: IT_UNTIL (u, x, xO , b) := ITERATE (IF (NOT b, u, x) , x, xO) 
28: IT_UNTIL (2 x, x, 3, 2 x > 100) 

29: 96 

30: ITS_UNTIL (u, x, xO , b) := ITERATES (IF (NOT b, u, x) , x, xO) 
31: ITS_UNTIL (2 x, x, 3, 2 x > 1000) 

32: [3, 6, 12, 24, 48, 96, 192, 384, 768, 768] 

6. Vielleicht haben sich andere DERIVE-User bereits mit diophantischen 

Gleichungen beschaftigt und konnen Herrn Lauer ihre Erfahrungen mittei- 
len? 



0 . 002 *c 


\ 


c c 

COS(t + 0.5) *SIN (2 • t ) , 0.005-c + SIN(t) -SIN(2-t) 

20 30 


c, 0, 84, 4 
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Berechnung und graphische Darstellung 
der Cornu-Spirale mir DERIVE 


Dr. F. Schumm, Stuttgart 

Bemerkung: Uber den Zusammenhang der Cornu' schen Spirale mit den Fresnel- 
schen Beugungserscheinungen, siehe Joos : Lehrbuch der theoretischen Physik, 
Seite 365 (1959) . (Diese Kurve fuhrt auch den Namen Klothoide. Sie ist eine 

transzendente ebene Kurve, deren Krummung proportional zur Bogenlange zu- 
nimmt und wird im Strahenbau verwendet da sie eine stetige Krummungsande- 
rung ermoglicht; siehe Wunderlich: Darstellende Geometrie 2, (1967)) 

Die Cornu' sche Spirale hat die folgende Parameterdarstellung mit dem Para- 
meter a: 


a a 

■ / ^ \ r 





' 2 N 






2 ) 




s • n 






s • n 

X (a) 


COS 


ds 

2 : 

Y (a) 


SIN 







- 

0 

2 J 




- 

0 

v 2 , 


Mit dem VECTOR-Kommando laht sich rasch eine Wertetafel ermitteln. Urn einen 
schnellen Uberblick zu erhalten wahlen wir die Grenzen 0 und 1 fur den Pa- 
rameter a mit einer Schrittweite 0.2: 

3: VECTOR ([a, X (a), Y (a)], a, 0, 1, 0.2) 

nach Simplify und geeigneter Einstellung von Precision (Approximate) und 
Notation (Decimal) erhalt man die folgende Tabelle: 








0 

0 

0 



0 

0 

0 


0.2 

0 .199 

-0 . 00418 



0.2 

0.199 

-0.00418 


0 . 4 

0.397 

-0 . 0333 



0.4 

0.397 

-0.0333 


0 . 6 

0 . 581 

-0 . 11 


TABLE([X(a), V(a>], a, 0, 1, 0.2) = 

0.6 

0.581 

-0.110 


0 . 8 

0 . 722 

-0.249 







1 

0 . 779 

-0 . 438 



0.8 

0.722 

-0.249 






1 

0.779 

-0.438 



Use the TABLE -function in DERIVE 5 


Urn die Punkte mit Plot direkt zeichnen zu lassen, darf die Tabelle nur die 
Werte von x und y, und nicht zusatzlich die Parameterwerte von a enthalten, 
sonst erhalt man die Fehlermeldung : "zuviele Variable". Daher loschen wir 

den Buchstaben a im Zeilenvektor [3:] . 

5: VECTOR ([X (a), Y (a)], a, 0, 1, 0.2) 

und erhalten nach neuerlichem Simplify: 



0 


0 

0 

.199 

-0 . 

00418 

0 

.397 

-0 

. 0333 

0 

.581 


0 . 11 

0 

. 722 

0 

.249 

0 

. 779 

-0 

. 438 


(TABLE( [X(a>, ¥{a)], a, 0, 1, 0.2)) COL [2, 3] 
0 0 
0.199 -0.00418 

0.397 -0.0333 


Again in Derive 5 


Diese Punkte konnen unmittelbar mit Plot ausgedruckt werden. Da dies noch 
zuwenige Punkte sind, muh man Schrittweite und Grenzen entsprechend anpas- 
sen. Bei der folgenden Kombination dauert die Berechnung allerdings auf ei- 
nem NEAT-AT einige Minuten : 








p 


8 


F. 
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7: VECTOR ( [X (a), Y (a)], a, -3, 3, 0.05) 

8: [[-0.605, 0.496], [-0.594, 0.448], [-0.562, 0.410], [-0.517, 0.390 

], [-0.467, 0.391], [-0.423, 0.413], [-0.392, 0.452], [-0.38~ 

~0 , 0 . 501] , [-0 . 388 , 0 . 549] , 


[0.388, -0.549], [0.380, -0.501], 

[0.392, -0.452], [0.423, -0.413], [0.467, -0.391], [0.517, -0.390] 

, [0.562, -0.410], [0.594, -0.448], [0.605, -0.496]] 


Das zugehorige Plot-Diagramm : 


y 




COMMAND : Center Delete Help Moue Options Plot Quit Range Scale Transfer 

Uindou aXes Zoom 


Plot in Derive 5 


Man kann naturlich auch direkt die Parameterdarstellung in folgender Weise 
eingeben und die Kurve dann mit Plot ausdrucken. Der Zeitbedarf ist aller- 
dings recht groft, bis die Graphik entstanden ist: 

9: [X (a), Y (a)] 




In spateren Derive - Versionen muss man die Datenpunkte nicht erst ausrechnen sondern wechselt 
gleich ins 2D-Plot-Fenster. Auflerdem konnten friiher die Punktgroflen nicht verandert werden. 
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Finanzmathematik mit 
Derive , Part 2 

Josef Bohm, Wiirmla 

3. Schuldtilgung (Tilgungsplane = Mortgage Tables) 

Wenn langerf ristige Kredite aufgenommen werden, ist es iiblich, sogenannte 
"Tilgungsplane" zu erstellen, aus denen der Stand der Schuldriickzahlung er- 
sichtlich ist. Fur das Ende jeder Riickzahlungsperiode (i.a. Jahre) werden 
ausgewiesen : 

- die aushaftende Restschuld (nach Zahlung der Annuitat) , [remaining 
debt ] 

- der Zinsenanteil , [payment of interest ] 

- der effektive Schuldendienst , - die Tilgungsquote, [redemption rate] 

- und die Annuitat. [ annuity ] 

Nach dem Modus der Riickzahlung unterscheidet man nun die Arten der Til- 
gungsplane. Hier sollen einige vorgestellt und mit DERIVE berechnet werden. 

3 . 1 Der unregelmafiige Tilgungsplan 


15. Eine Schuld von 150000.- wird aufgenommen. Als Riickzahlungen werden 
fur die ersten beiden Jahre 20000.- und 40000.- und fur das 4. Jahr 
65000.- vereinbart. Mit dem 6. Jahr soil die Riickzahlung abgeschlossen 
sein . 

Wie lautet der Tilgungsplan bei 10% Zinsen? 

A debt for 150000.- is borrowed. As repayments have been arranged 
20000.- and 40000.- for the first two years and 65000.- for the 4th 
year. The paying back is to be settled within six years. What's the 
redemption table at an interest rate of 10%? 


Wie bei der Investitionsrechnung definieren wir einen Vektor, in dem der 
Reihe nach zuerst das Zahlungsjahr und dann die Annuitaten angegeben wer- 
den. ( payments and when they are due are given in a vector) 

1 : paym : = 

DIM (paym) Ab DERIVE 2.xx werden konstante 

2 : n : = + 1 

2 GroBen mit Kleinbuchstaben bezeich- 

3 : r (i) := 1 + i% ne. 

4: ann(j) := ELEMENT (paym , 2-j) 

5: time ( j ) := ELEMENT (paym , 2-j - 1) 

Mit einfachen Formeln der Finanzmathematik wird nun die Schluhannuitat , das 
ist der Zahlungsrest am Ende der Laufzeit berechnet: 

n n - 1 n - time ( j ) 

6: last_ann(k, i) := k-r(i) - Z ann(j)-r(i) 

j=l 

REM_DEBT ( k, i , h) ist die Restschuld nach Bezahlung der Annuitat im h-ten 
Jahr der Riickzahlung: (last annuity and remaining debts) 

h h h - j 

7: rem_debt(k, i, h) :=k-r(i) - Z ann(j)-r(i) 

j=l 

8: ints(k, i, h) := rem_debt(k, i, h - 1 ) • i% 
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In einer Zeile des Tilgungsplanes werden der Reihe nach ausgewiesen: 


Jahr, 

Restschuld, 

Zinsendienst , 

Schuldendienst , 

Annuitat 

Year, 

Balance, 

Interests , 

Principal , 

Annui ty 

(h) 

Rh) 

(Zh) 

(th) 

(Ah) 


Anschliehend wird wieder der so universell zu verwendende VECTOR-Bef ehl 
eingesetzt, um den Tilgungsplan - aber nur bis zur vorletzten Zeile - zu- 
sammenzuset zen . 


9: row_(k,i,h) : = [h , rem_debt ( k , i , h ) ,ints(k,i,h) , ann (h) -ints ( k , i , h) ,ann(h) ] 

10: titel := [ ["h" , "Rh" , "Zh" , "th" , "Ah"] ] 

11: red_tab ( k , i, m) := APPEND (titel, VECTOR ( row_ ( k , i, h) , h, m) ) 


(In Derive 2.xx hiefi es noch ROW, aber in Version 5 ist ROW eine Funktion, APPEND musste damals als 
APPEND VECTORS aus einer Utility-Datei geladen werden.) 

Die Terme 1: bis 11: konnten in einer Datei RED_TAB.MTH Oder TPL.MTH abge- 
legt werden. 

Nun zu Beispiel 15: 

12: paym := [l, 20000, 2, 40000, 3, 0, 4, 65000, 5, 0] 

13: last_ann (150000 , 10) 

14: 96309.95 

Wir erweitern die Liste paym umd das Jahr 6 und die Schluhannuitat und kon- 
nen mit red_tab (Kapital , Zinsfu(5 , Lauf zeit ) den kompletten Tilgungsplan abru- 
f en . 

15: paym := [l , 20000, 2, 40000, 3, 0, 4, 65000, 5, 0, 6, 96309.95] 

16: red_tab (150000 , 10, 6) 


h 

Rh 

Zh 

th 

Ah 

1 

145000 

15000 

5000 

20000 

2 

119500 

14500 

25500 

40000 

3 

131450 

11950 

-11950 

0 

4 

79595 

13145 

51855 

65000 

5 

87554.5 

7959.5 

-7959.5 

0 

6 

0 

8755.45 

87554.5 

96309.95 


Diese Behandlung hat noch zwei Schonheitsf ehler : 

die besondere Berechnung der Schluhannuitat und ihre nachtragliche Einbin- 
dung in den Zahlungsvektor sowie die Nichtbeachtung einer allfalligen Run- 
dung . 

Die neue Version von DERIVE macht nun beides moglich: 

Die Rundung wird im letzten Beispiel dieses Beitrags demonstriert . Der Run- 
dungsalgorithmus auf n Nachkommastellen funktioniert f olgendermahen : 

7/ x INT(x*10” +0,5) 

rci (x) = 

10 ” 

Dabei bedeutet INT(a) die grohte ganze Zahl < a. Diese Funktion wird von 
DERIVE in der Datei MISC.MTH unter dem Namen FLOOR (x,l) gemeinsam mit 
MOD (x , 1 ) bereit gestellt. ( Jetzt sind FLOOR und MOD bereits implementiert. ) 
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In der Utility-Datei VECTOR. MTH findet man nun auch die Funktion 
APPEND_VECTORS (u,v) r mit der man Matrizen und Vektoren zusammenhangen kann . 
(Das 1st historisch , jetzt wird mit APPEND und ohne Nachladen gearbeitet.) 

Wir laden diesen Befehl aus VECTOR. MTH mit Merge zu den Ausdrucken von vor- 
hin. Auherdem habe ich die mdglicherweise nutzliche Funktion: 

PART_RT (Kapital , Zinsfup) , Anf angszeile , Endzeile ) 

definiert, mit der man den Tilgungsplanauch teilweise ausgeben kann. 


18 : 

part_rt ( 

k, 

i, u , 

v) := VECTOR ( row_ ( k , i, h) , 

h, 

u , v) 

19: 

red_tab ( 

k, 

i) : = 

APPEND ( titel , part_rt(k, i, 

1, 

n - 1) , 


[[n, 0, 

rem 

_debt ( 

k, i, n - 1) , rem debt(k, i 

, n 

- 1) -i% 


rem_debt 

(k, 

i, n 

- l).r(i)]]) 



20 : 

paym : = 

[1, 

20000 

, 2, 40000, 3, 0, 4, 65000, 

5 , 

0] 


Die Parameter k und i lassen sich auch als Konstante definieren. 


21 : 

[k : 

= 150000, i := 

10] 

22 : 

row_ 

( k , i , 3) 


23: 

[3, 

131450, 11950, 

-11950 , 

24 : 

part 

_rt ( k , i, 2, 5) 



2 119500 14500 25500 40000 

3 131450 11950 -11950 0 

25: 4 79595 13145 51855 65000 

5 87554.5 7959.5 -7959.5 0 


Nun sofort und ohne weiteren Zwischenschritt zum kompletten Tilgungsplan : 
26: red_tab(k, i) 



Und so andert sich der Tilgungsplan, wenn der Zinsfuh urn zwei Prozent 
steigt, Oder urn ein Prozent fallt: 

27: red_tab ( k , 12) 

28: red_tab(k,9) 


h 

Fth 

Zh 

th 

Ah 


h 

Fth 

Zh 

th 

Ah 

1 

148000 

18000 

2000 

20000 


1 

143500 

13500 

6500 

20000 

2 

125760 

17760 

22240 

40000 


2 

116415 

12915 

27085 

40000 

3 

140851.2 

15091.2 

-15091.2 

0 


3 

126892.3 

10477.35 

-10477.35 

0 

4 

92753.34 

16902.14 

48097.85 

65000 


4 

73312.66 

11420.31 

53579.68 

65000 

5 

103883.7 

11130.40 

-11130.40 

0 


5 

79910.80 

6598.139 

-6598.139 

0 

6 

0 

103883.7 

12466.04 

116349.7 


. 6 

0 

79910.80 

7191.972 

87102.77 
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3 . 2 Die Annuitatenschuld 


16. 85000.- sind bei 8,75% Zinsen innerhalb von 10 Jahren durch gleiche 

Raten zuruckzuzahlen . 


a) Wie lautet die 7. Zeile des Tilgungsplans? 

b) Erstelle die 3. - 7. Zeile des Tilgungsplans. 

c) Wie lautet der komplette Plan? 

d) Wie andert sich die 7. Zeile, wenn der Zinsfut um 1% erhoht wird? 

85000.- are to be paid back within 10 years at 8.75 % interest by equal 
annul ties . 

a) What is the 7th row of the redemption table? 

b) Give the rows 3 to 7 of the redemption table! 

c) Give the complete table! 

d) How does row 7 change , if interest is increased by 1%? 


Die Terme bis 7: konnten (ohne Kommentarzeilen) in der Hilfsdatei ANN.MTH 
zusammengef asst werden. Aufgabe 16 wird mit dieser ANN . MTH-Datei gelost: 

Annuitatenschuld: Kapital (Amount) K, Zinsfuft (Interests) i and Zeit (Years) n. 

#1: r(i) := 1 + ±y. 

Berechung der Annuitat (Annuity) 

hi/. 

a(k, i, n) := 

#2: -n 

1 - r(i> 

oder -pmt(i%,n,k) 

Berechnung der Tilgungsquoten fur das Jahr h (Principals for year h) 

h - 1 

#3: i, n, b) := (a(k, i, n) - k-ix)-r(i) 

Restschuld zu Beginn des h-ten Jahres 
Balance at begin of year h 

h 

r_r(k, i, n, l)-(r(i) - 1) 

#4: ren_debt (k, i, n, b) := k - 

i /. 

oder k + FVAL(i% ,h ,r_r(k ,i ,n ,1 )) 

Ausgabe der h-ten Zeile des Tilgungsplans 

#5: rou_(k, i, n, b) := [b, rem_debt (k, i, n, h), a(k, i, n) - r_r(k, i, n, h), r_r(k, 

i, n, h), a{k, i, n)] 

#6: titel : = [[b, Rb, Zb, th, fib]] 

Ausgabe des Plans von Zeile u bis Zeile v 
Output of the plan from row #u through row#v 

#7: part_rt{k, i, n, u, u) := flPPEND{t itel, UECTOR{row_{k, i, n, b), b, u, u)) 

#8: Notation := Decimal 

#9: Notat ionDigits := 7 

Kompletter Plan 

#10: red_tab{k, i, n) := part_rt{k, i, n, 1, n) 

#11: [i := 8.75, k := 85000, n := 10] 

#12: rou_(k, i, n, 7) 

[7, 33306.53, 3733.804, 9365.516, 13099.32] 


#13: 


D-N-L #2 


J. Bohm: Financial Mathematics 2 


P 13 


#14: 

part_rt(k, i. 

n, 3, 7) 





h 

Fth 

Zh 

th 

Ah 



3 

66484.96 

6403.333 

6695.987 

13099.32 



4 

59203.07 

5817.434 

7281.886 

13099.32 

#15: 


5 

51284.02 

5180.269 

7919.051 

13099.32 



6 

42672.05 

4487.351 

8611.968 

13099.32 



. 7 

33306.53 

3733.804 

9365.516 

13099.32 

#16: 

red_ 

tab(k, i. 

n) 





h 

Fth 

Zh 

th 

Ah 


#17: 


1 79338.17 7437.5 

2 73180.94 6942.090 

3 66484.96 6403.333 

4 59203.07 5817.434 

5 51284.02 5180.269 

6 42672.05 4487.351 

7 33306.53 3733.804 

8 23121.53 2914.322 

9 12045.35 2023.134 

10 0 1053.968 


#18: rou_(k, 0.75, n, 7) 




5661.820 13099.32 

6157.230 13099.32 

6695.987 13099.32 

7281.886 13099.32 

7919.051 13099.32 

8611.968 13099.32 

9365.516 13099.32 

10184.99 13099.32 

11076.18 13099.32 

12045.35 13099.32 

26170.13, 260.7296, 8593.824, 8854.554] 


Die Funktion PART_RT ( k, i , n, u , v) ist dann praktisch, wenn der Tp . seines Um- 
fangs wegen am Bildschirm und damit auch bei der Druckausgabe nicht in die- 
ser Matrixform ausgegeben wird. Mit PART_RT kann man ihn teilweise ermit- 
teln und stuckeln. 


3 . 3 Die Prozentschuld 

Hier ist neben dem Kapital K und der 
Verzinsung i die Hohe der jahrlichen 
konstanten Ruckzahlungen , - der An- 
nuitaten - , gegeben . Vorrangig ist 
die Laufzeit zu ermitteln. 

Der Name „Prozentschuld" ruhrt daher, 
dass die Hohe der Annuitat oft als 
Prozentanteil der Gesamtschuld defi- 
niert wird. 

Ich habe diesen Tilgungsplan auf drei 
Arten erstellt: 

Zuerst fur DERIVE-User, die - noch - 
nicht im Besitz der neuesten Version 
sind. Hier muss man ahnlich wie beim 
allgemeinen Tilgungsplan in 3.1 nach 
Berechnung der Laufzeit die letzte 
Zeile separat ermitteln. 



VECTOR ( [ (3 . 5 + c • SIN ( 4 • t) ) COS (t) , 

(3 . 5+c • SIN ( 4t ) ) SIN (t ) ] ,c, -2,2,0.25) 


Dann zeige ich mit Verwendung der INTEGER-Funktion und APPEND_VECTORS eine 
wesentlich komf ortablere Losung und schliehlich will ich auch einen kleinen 
Vorgeschmack auf die Programmierbarkeit von Derive 2 . xx geben . Ja, Sie ha- 
ben richtig gelesen: Derive 2.xx ist programmierbar ! Auherdem werde ich in 
der letzten Variante den schon besprochenen Rundungsalgorithmus einbauen: 
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17. Die Schuld von 120000.- ist bei i=10% durch Annuitaten von 15600.- zu 
tilgen . 


a) Wie lautet die 10. Zeile des Tilgungsplans? 

b) Bestimme die Gesamtlauf zeit . 

c) Wie heiftt die letzte Zeile? 

d) Erstelle den Tilgungsplan von der 8. bis zur 13. Zeile. 

A debt of 120000 . - at 10.5 % interest is to be paid back by annuities 
of 15600.-. 

a) What is the 10th row of the redemption table? 

b) Calculate the running time of the pay-off . 

c) Give the last row. 

d) Give the schedule of rows 8 through 13. ? 


Hier ist die Datei PROZSCH.MTH dokumentiert . Der Formelbereich geht bis zum 
Ausdruck # . . : 


Prozentschuld : Kapital k, Annuitat a, Zinsfuh i (in %) 


#1: r ( i ) : = 1 + i% 


LN 


k • i° 


1 - 


#2: n(k, a, i) : = 1 + 


-1 

LN ( r ( i ) ) 


Berechnung der Tilungsquoten und Reschschulden 
Calculation of principals and the balances 

h - 1 

#3: r_r(k, a, i, h) := (a - k«i%)«r(i) 


#4: rem_debt(k, a, i, h) 


r r ( k , a , i , 1 ) 


k - 


i) 


Ausgabe der Zeile h (nicht der letzten) 

Output of row #h (not the last one) 

#5: row_ ( k , a , i , h) := [h , rem_debt ( k , a , i , h ) , rem_debt ( k , a , i , h - 1 ) • i% , 

r_r ( k , a , i , h ) , a] 

Ausgabe letzten Zeile (m = INT (n) ! ) 

Output of last row (m = INT (n) ! ) 

#6: lrow(k,a,i) := [m , 0 , rem_debt ( k , a , i ,m - 1 ) • i% , 

rem_debt ( k , a , i , m - 1), rem_debt ( k , a , i , m - 1 ) • r ( i ) ] 

Ausgabe der Zeile h (nicht der letzten) 

Output of row #h (not the last one) 

#7: titel := t ["h","Rh","Zh","th","Ah"]] 

#8: part_rt ( k , a , i , u , v) := APPEND ( titel , VECTOR ( row_ ( k , a , i , h ), h , u , v) ) 

Kompletter Plan ohne letzte Zeile / complete table without last row 


#9: 

red_tab ( k , a , i ) 

I := APPEND 

(titel 

, part_rt ( k , a , i , 1 , m- 1 

#10: 

[k := 120000, 

H- 

II 

M 

O 

Cn 

a : = 

15600] 

#11: 

row_ ( k , a , i , 

10) 



#12 : 

[10, 71026. 26 : 

, 8231.454, 

7368 . 

545, 15600] 
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#13: 

m = 

n ( k , a , i ) 




#14: 

m = 

17 . 51212 




#15: 

lrow 

(k, a, i, 17 

) 



#1 6 : 

[17, 

0, 777.6281, 7405.982, 8183.610] 

#17: 

part_ 

-H 

-p 

p 

i 

8, 13) 





h 

Rh 

Zh 

th 

Ah 



8 

85063 . 17 

9565.278 

6034 . 721 

15600 



9 

78394 . 80 

8931 . 633 

6668.366 

15600 

#18: 


10 

71026.26 

8231 . 454 

7368 . 545 

15600 



11 

62884 . 01 

7457 . 757 

8142.242 

15600 



12 

53886.84 

6602 . 821 

8997 . 178 

15600 



13 

43944 . 95 

5658 . 118 

9941 . 881 

15600 


Nun die aktuellere Variante: 

Die Zeilen 1: bis 4: sind voranzustellen; damit andern sich auch die 
Ermittlung von Laufzeit und letzter Zeile. Eine "Nullzeile" wird 
eingefuhrt, um Tilgungsplane uber ihre tatsachliche Laufzeit ausge- 
ben zu konnen, ohne vorher die Laufzeit berechnen zu miissen. 


18. 50000 sind durch Annuitaten von je 15000 bei 9% Zinsen zu tilgen. 

Wie lauten die Zeilen 23 bis 30 des Tilgungsplans ? 

Wie lange dauert die Riickzahlung? 

150000.- must be paid back by annuities of 15000.- at 9 percent inter- 
rests . 

What are the rows 23 to 30 of the payback table? 

What is the running time of the pay-off? 


l: 


b ATAN [COT [ * b a ]] 


MOD (a, b) := -| 


2: 

FLOOR (a, b) := a 

MOD (a, b) 
b 

3: 

INTG (x) := FLOOR (x, 

1) 

4: 

APPEND_VECTORS (u, v) 

VECTOR (IF (ra_ < DIMENSION (u), ELEHEN' 


~T (u, m_), ELEMENT (v 

, m_ - DIMENSION (u))), m_, DIMENSION (u) ^ 


DIMENSION (v)) 


5: 

R (i) := 1 + i% 



This is history, all functions are full 
implemented now. It is not necessary 
to preload them from any utility files. 

APPEND_ VEC T ORS became 
APPEND. 


ttl : r(i) := 1 + iY. 

b - 1 

#2: r_r(k, a, i, b) : = (a - k-iK)-r(i) 

#3: ren_debt{k, a, i, b) := k - FUAL{ix, Ji, - r_r{k, a, i, 1)) 

#4: rou_(k, a, i, b) := [h, rem_debt (k, a, i, b), a - r_r(k, a, i, b), r_r(k„ a, i, b), 

a] 

#5= Berecbnung der Laufzeit: 

#6: n(k, a, i) := FLOOFt( NPER< i'A , -a, k)) + 1 

#7: Ausgabe der letzten Zeile / last row 

#8: lrow{k, a, i) := [n{k, a, i), 0, ren_debt{k, a, i, n{k, a, i) - l)-itf, ren_debt{k, 

a, i, n(k, a, i) - 1), rem_debt (k, a, i, n(k, a, i) - l)-r(i)] 

#9: Nullzeile / Zero line 

z(b> := [b, 0, 0, 0, 0] 


# 10 : 
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#12: kompletter Plan / full table 

#13: red_tab(k, a, i) := APPEND(t itel, UECTOR(row_(k, a, i, b), b, 1, n(k, a, i) - 1), 

[lrow(k, a, i>]) 

#14: erweiterter Plan bis sur Zeile u > n / extended table up to row ttu > n 

#15: rtx(k, a, i, u) := APPEND(red_tab(k, a, i), UECTOR<s<b), b, n(k, a, i) + 1, u)) 

#16: Teilt ilgungsplan Zeile ttu bis ttu / Partial table row ttu to ttu 

#17: part_rt(k, a, i, u, u) := APPEND(t itel, U ECT OR (U ECTOR (ELEMENT (rtx(k, a, i, u), j, 

b), b, 1, 5), j, u + 1, u + 1>) 

#18: part _rt (150000, 15000, 9, 23, 30) 


b 

Rb 

Zb 

th 

Ah 


23 

45702.09 

5012.099 

9987.900 

15000 

Die Ruckzahlung dauert 

of f ensichtlich 27 Jahre . 

24 

34815.28 

4113.188 

10886.81 

15000 

It needs 21 years to pay 

25 

22948.65 

3133.375 

11866.62 

15000 

back . 

26 

10014.03 

2065.379 

12934.62 

15000 


27 

0 

901.2631 

10014.03 

10915.29 

Below is part rt in the modern 

28 

0 

0 

0 

0 

version. 

29 

0 

0 

0 

0 


30 

0 

0 

0 




part_rt(k, a, i, u, u) := APPEND(t itel, UECTOR((rtx(k, a, i, u)) , j, u + 1, u + 1)) 
#20: j 

Wenn wir auf 2 Dez.stellen runden wollen, erganzen wir den Funkti- 
onsvorrat mit der Funktion RD(x) . ROW_ ( k, a, i , h) und LROW(k,a,i) wer- 
den entsprechend adaptiert. APPEND_VECTORS kann nun weggelassen wer- 
den . An seiner Stelle wird ja n programmiert " . 

FLOOR <100 x +0.5) 

#21: rd(x) := 

100 

#22: row_(k, a, i, b) := [b, rd(ren_debt (k, a, i, b)), rd(a - r_r(k, a, i, b)), 

rd(r_r(k, a, i, b)), rd(a)] 

#23: lrou(k, a, i) := [n(k, a, i), 0, rd(ren_debt (k, a, i, n(k, a, i) - i)-ix), 

rd(ren_debt (k, a, i, n(k, a, i) - 1)), rd(ren_debt (k, a, i, n(k, a, i) - l)-r(i))] 

Hier wird mit der in vielen ho- 
heren Programmiersprachen ubli- 
chen IF-Abfrage entschieden, ob 
noch eine Normalzeile ROW_, die 
letzte Zeile LROW oder eine 
uberflussige Nullzeile Z aus- 
zugeben ist. 

Like in many programming lan- 
guages we use an IF-construct 
to decide whether a standard 
ROW_, the last row LROW or a 
zero row Z must be produced . 


VECTOR ( [ (3 . 5+cSIN (4t) ) COS (-t) , (3 . 5+cSIN (4t) ) SIN (2t) ] f 
c, -2,2, 0.5) 



Compare with graph on page 13! 
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Lassen Sie diese Konstruktion auf Sie einwirken. Wir werden sicher 
in den nachsten Ausgaben des D-N-L eingehend auf die neuen Moglich- 
keiten von Derive 2 . xx eingehen. 


#24: part_rt(k, a, i, u, y) := APPEND(t itel, UECTOR(IF(Ji < n(k, a, i), row_(k, a, i, Ji), 

IF(Ji = n(k, a, i), lrou(k, a, i), 2 <h))), h, u, y)) 


#25: 

red_tab(k, a, i) := part_rt(k, a, i 

, 1, n(k. 

■rl 


#26: 

part _rt {150000, 15000, 9, 

20, 27) 






h 

Rh 

Zh 

th 

flh 



20 

73259.82 

7287.51 

7712.49 

15000 



21 

64853.2 

6593.38 

8406.62 

15000 



22 

55689.99 

5836.79 

9163.21 

15000 

#27: 


23 

45702.09 

5012.1 

9987.9 

15000 



24 

34815.28 

4113.19 

10886.81 

15000 



25 

22948.66 

3133.38 

11866.62 

15000 



26 

10014.03 

2065.38 

12934.62 

15000 



27 

0 

901.26 

10014.03 

10915.3 


AbschlieLend sei noch gesagt, daft dies naturlich nur Vorschlage fur 
den Einsatz von DERIVE sein konnen. Da gibt es sicher elegantere 
Moglichkeiten und Erganzungen. Ich wollte mit meinem Beitrag nur 
zeigen, daft auch f inanzmathematische Anwendungen mit DERIVE rasch 
und erfolgreich behandelt, und daft auch sehr mathematischen Bezuge - 
hier zu Vektor- und Matrizenrechnung - hergestellt werden konnen. 

Ich freue mich auf Ihre diesbezuglichen Erf ahrungsberichte oder zu- 
satzlichen Beitrage . 
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Ubernahme von Texten und Graphiken aus 
DERIVE in Word5 mit Hilfe von 
CAPTURE.COM 


Dr. Felix Schumm, Stuttgart 

1. Vorschlag zur Installation von DERIVE und starten iiber eine BATCH-Datei. 

Im Heft Nr.l der D-N-L fragt Herr Dr. G. Bragard, Aachen an, wie man DERIVE mit einer Hercules- 
karte so installieren kann, dafi man nicht bei jedem Neustart HGC an-wahlen muB. Hier ein Losungs- 
vorschlag. Die anschlieBenden Hinweise zur Ubernahme von Texten und Graphiken in W0RD5 be- 
ziehen sich auf diese Installation. (Sie kann leicht auf spezielle Bediirfnisse angepaBt werden. Anmer- 
kung von D-N-L) 

Ich ordne auf der Festplatte C: Verzeichnisse und Dateien in folgender Weise an: 

Stammverz. Unterverz. Dateien 

root subdirectories files 

Vc:\ 

config.sys 

autoexec.bat 


\BATCH\ 


DERIVE . BAT 
WORD . BAT 

{weitere . BAT-Dateien } 


\DERIVE\ 


DERIVE . EXE 

DERIVE . INI 

HERCULES.COM 

{weitere DERIVE-Dateien } 


\WORD\ 


WORD . EXE 

CAPTURE . COM 

CAPTURE . INI 

{weitere WORD-Dateien } 


\DOKU\ 


TEXT . TXT 
BILD. SCR 

{weitere Texte und Bilder fur WORD} 


Aufbau der Batch-Dateien: 


Zuerst wird der Suchpfad in der AUTOEXEC.BAT entsprechend adaptiert: 
PATH C:\DOS; C:\BATCH; C:\WORD;... 
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WORD . BAT 

DERIVE . BAT 

echo off 

echo off 

els 

els 

c : 

c : 

cd c:\doku 

capture / a=0 6/ t=N/r=N/i=J/b=J 

word/h 

cd derive 

c : 

hercules 

cd\ 

derive 

els 

c : 


cd\ 

els 


Wenn man keine HERCULES Karte, sondem einen CGA-, EGA- oder VGA-Schirm verwendet, fallt 
das /h in WORD.BAT, sowie hercules in DERIVE.BAT weg. Informieren Sie sich bitte im WORD- 
Handbuch iiber CAPTURE.COM! 

Nach Eingabe von ’’derive” erscheint zun„chst das Menu des Programms CAPTURE. Man andert oder 
bestatigt die Vorgaben, indem man die Q-Taste driickt. Danach erscheint das DERIVE-Menii. Ich 
arbeite grundsatzlich im Graphik-Modus. Beim ersten Start wahle man unter 

Options Display die Moglichkeiten: graphics, high, hercules (oder entsprechend) und extended. 

und unter Transfer Print Options die Einstellungen: all, extended und height=50. 

Dann speichert man diese Einstellungen mit Transfer Save State in die Datei DERIVE.INI ab. Bei 
erneutem Start von DERIVE ist diese Vorwahl dann automatisch eingestellt. 

( Anmerkung von D-N-L: Es lassen sich verschiedene .INI-Dateien fur unterschiedliche Einsatzmog- 
lichkeiten speichern und mit Transfer State Load aufrufen (z.B.: MONO.INI, INVERS.INI, ...) 

2. Ubernahme des Algebra-Fensters von Derive in Word5. 

Dies ist unkritisch. Nach der Befehlsfolge Transfer Print File gibt man den Dateinamen, sowie den 
Pfad ein, in dem WORD den Text suchen soil (z.B.: C:\DOKU\NAME.TXT). Unbedingt sollte man 
jedoch, wie zuerst beschrieben die Einstellungen extended und height=50 eingetragen haben, um ei- 
nerseits die Sonderzeichen des ASCII-Satzes, die die meisten Drucker unterstiitzen zu erhalten, und 
um andererseits groae Integrale oder auch Doppelbriiche und Summen in Originaldarstellung aufs 
Papier zu bringen (und nicht als INT(. ...)). Beim Bearbeiten der Texte mit WORD5 sollte man minde- 
stens bei Briichen den Zeilenabstand auf 0,5 setzen (FORMAT-ABSATZ) und mit FORMAT- 
ZEICHEN eine kleinere SchriftgroBe wahlen. 

(Es bewahrt sich auch, den Bereichsrand so weit wie moglich nach auBen zu setzen, und die ersten 
Spalten mit SHIFT F6 zu entfernen.) 

3. Ubernahme des PLOT-2D/3D Fensters von Derive in Word5 mit Hilfe von CAPTURE.COM 

Wenn man SHIFT-DRUCK (PrtSc) driickt, erscheint links oben am Schirm die Meldung Dateiname: 
captOOl .scr. Hier geben wir den Namen und Pfad ein, aus dem WORD die Datei wieder finden soil 
(z.B.: C:\DOKU\KURVEl.SCR). Nach ENTER erscheint ein Rahmen, den man mit den Pfeil- und der 
Tabulatortaste verschieben kann, der jenen Bereich anzeigt, derin ein fur WORD lesbares Bild- 
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format iibertragen werden soil. 1st der Ausschnitt passend gewahlt, bestatigt man neuerlich mit 
ENTER und der Vorgang ist abgeschlossen. 


In WORD5 laBt sich nun an der Stelle, an der der Cursor gerade steht, die Graphik mit der Befehlsfol- 
ge BIBLIOTHEK-VERKNUPFEN-GRAPHIK einbinden. Man muB nur noch den Namen der Bildda- 
tei (z.B.: KURVE1.SCR) eingeben. Gegebenenfalls kann man die GroBe der Graphik noch anpassen 
und diese mit FORMAT-RAHMEN (Rahmen)noch umranden. 


Mit der Layoutkontrolle (Ctrl+F9) laBt sich die Plazierung noch kontrollieren. 

Die Bilder zur Cornu-Spirale in dieser Ausgabe des D-N-L sind auf diese Weise nach WORD iibertra- 
gen worden. 


4. Additional Comments from J. Bohm 

There are some other possibilities to print, to save and/or to manipulate DERIVE graphics. PIZZ- 
AZZ-PLUS is a fine program to create first class hardcopies. Many options make it easy to change, to 
rotate, to smooth, to shrink and to save the generated graphics. 

If I want to change a DERIVE-plot - in order to add a description and/or a scale, I prefer working with 
FRIEZE, a PC-Paintbrush utility. Once installed FRIEZE remains memory resident until I re -boot the 
computer or turn off the power. FRIEZE allows to capture and save graphics images from other pro- 
grams in PC-Paintbrush ? s PCX-format. WORD5 is able to read this format, as it does with its own 
.SCR-files mentioned above. 

To show an example - having in mind the Letter of the Editor - 1 captured the following family of 
curves first by using CAPTURE.COM and transferred it to a WORD'S .SCR-file. Then I used FRIEZE 
to save the same graphics for later "embellishing" with PC-Paintbrush. Please admire the result of my 
efforts below! Thank you so much! 

All the graphics in this issue of the D-N-L were produced in one of these ways. The time of cut and 
paste has gone by. 


VECTOR ([ (3. 5+c SIN ( 4 t ) ) COS ( -t ) , (3.5+c SIN(2t)) SIN (2t) ] , c, -4, 0) 



That's the original, captured by 
WORD'S CAPTURE.COM. 

(This is history!!) 
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Felix, the Cat 


And this is the 
paining after 

working out the 
DERIVE-plot with 
PC-Painbrush. 
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Solving Odes Using 
DERIVE ' s ODE1.MTH 


Josef Bohm, Wurmla 


The file 0DE1.MTH is helpful to solve first-order ordinary differential 
equations. I want to explain the capabilities of this file, complete it and 
make its use a bit easier. 


The simplest Odes are those which can be solved by separating the vari- 
ables. Let's look at an example: 


2 

Example 1 : ( 1 + y ) dx + xy dy = 0 

Give the general solution and give the curve of solution con- 
taining then point P ( 3 / 2 ) . 

Sketch the direction field, given by (1) , some curves of the 
general solution and the curve of the special solution! 


d y 

It is necessary to bring the given Ode to the form: = p (x) 

dx 

to use the function SEPARABLE: 


1: SEPARABLE (p, q, x, y, xO, yO ) 


y 

yO 


1 


q 


X 


dy = J 

xO 


p dx 


After separating the variables I get: p 

and this leads to: (in Derive 5!) 


1 1 + y2 

and q = 

x y 


q(y) 


#2: Logarithm := Expand 


#3: SEPARABLE 


1 1 + y 

x y 


2 2 
LN(y + 1) LN(y0 + 1) 


= LN(xO) - LN(x) 


#4: [x :e Real (0, ®), y :e Real (0, «»)] 


#5 : SOLUE 


2 2 
LN(y + 1) LN(y0 + 1) 


LN(x0) - LN(x), y 


2 2 2 2 2 2 
\T(x0 (y0 + 1) - x ) \T(x0 (y0 + 1) - x ) 

#6 : y = - v y = 


tt7: SEPARABLE 


1 1 + y 


r x j. y, 3 j- 2 


LN(y + 1) LN{5) 


2 


= LN(3) - LN(x) 


tt8 : SOLUE 


2 


LN(y + 1) LN(5) 


= LN{3) - LN(x)„ y 


2 2 
T{45 - x ) T{45 - x ) 

tt9 : y = - v y = 


2 x - 45 

ttl0: y = 


x 
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With feld() you can vary the length of the tangent pieces by setting the 
values for parameter t, -0 . 25 < t < +0 . 25 is recommended. Plot the direc- 
tion field in blue. 


Set the appropriate scale, switch off the automatical change of colors, 
switch on Approximating Before Plotting and plot .... and then you will 
watch how the direction field will spread across the defined plane. 

I've brought the function term in a more convenient form and build up two 
vectors of the family of curves, which represent the general solution. 


#24: 


U ECTOR 


2 2 
I(c - x ) 

, c, -20, 20, 2 

x 


#25: 


U ECTOR 



c, -20, 20, 2 


Plot them all in red 



and finally add the special solution with the given point: 


# 36 : 


2 

V 


2 

x - 45 

, [3, 2] 

2 

x 


Many things have changed since 1991: You need not loading the utility files. Derive recognizes all 
functions which are provided in the utility files (which you can find in the MATH-f older. It was not 
possible to plot implicitly given functions (#36). 


We can do this now, see expression #38 . 


2 

V = 


2 2 
c - x 


2 


#38: U ECTOR 


x 


-, c, —20, 20, 2 
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Example 2 : 

(1 + e x ) y y 1 = ex 

spec. sol. P(l;l) 

Example 3 : 

y T sin x = y In y 

y(n/2) = 1/2 


#26: spec_sol_sep 


1, 1 


1 + e 


x y 


#27: — — = LN(e X + 1} - LN(e + 1) 

2 2 


#28 : SOLUE 


2 

y 


1 x 

= LN(e + 1) - LN(e + 1}, y 

2 2 


X X 

#29: y = - sf(2 LN(e + 1) - 2 LN(e + 1) + 1) v y = vf(2 LN(e + 1) - 2 LN(e + 1) + 1) 

x 2 

#30: (y = - \T{2 LN(e + 1) - 2 LN(e + 1) + 1>> 


2 x 

#31: y = 2 LN(e + 1) - 2 LN(e + 1) + 1 


#32: f eld 


x, y, 0, 5, 0.5, -3, 3, 0.5 


V <e + 1) 



#34: spec_sol_sep 


1 n 1 

, y LN(y), , 

SIN(x) 2 2 


LN(y) 

#35: LN | 1 - n i = LN|TAN 

LN{2) 


■H) 


f))- 


f f LN<y> 1 

#36: SOLUE LN - n i = LN |TAN 

l l LN<2) J 

f f f LN<y> 

#37: SOLUE SOLUE LN 1 - n-i = LN |TAN 

l L L LN<2) 

#38: x : £ Complex 

#39: y : £ Complex 


f))- 


y I - y 


#40: y = IF |-n < IM |lN |lN( 2) -TAN |-^-J J + n-ij < n, If|-h < Im|- LN(2) TAn|-^-JJ < n, 2 

]] 


TAN{x/2)l 


#41: gen_sol_sep 


SIN(x) 


y-LN(y) = LN(LN(y)) - LN TAN 


t]] = 


#42: SOLUE |LN(LN<y>) - LN [tAN [-^] ] = c, yj 
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#43: y 


IF 


-n < IM LN TAN 


'll 


+ c 



f c f x V 

| e C TAN(x/2) 

n, IF 

-n < IM e TAN — 

^ n, e 


c 

#44: ^ TAN(x/2) 

y = e 


#45: 

#46: 

#47: 


SOLUE 


1 TAN(n/2/2) 

= e 

2 


c = LN(LN{2>) + 3 n i v c 


e TAN(x/2) 
SUBST ly = e , c. 


= LN(LN{2>) - n i v c = LN(LN(2)> + n i 

- TAN(x/2) 

LN(LN(2>) + n ij = (y = 2 ) 


- TAN(x/2) 

#48: y = 2 


It's interesting to compare De- 
rive 's solution from 1991 with its 
treatment in 2004. 

Domain declarations are very 
important. 


SPEC_SOL_SEP C-sIsSo"' y LN (y) ’ -T' -h 

111 LN (2) (COS^(x) - 1) 3 “ LN C SIN^x) 1 " * * = 0 
.LN (- LN (y) / (LN (2) (COS (x) - 1))) - LN (1 / SIN (x)) - i Z 
Ztt = 0 

SIN (x) LN (y) 

LN (2) (COS (x) - 1) " 1 


. y = 2 COT (X) - 1 / SIN (x) 


Example 4: The logistic growth curve 

Logistic growth is described by a difference equation: 
Increase of a population y during a time interval dx is 
proportional to its present quantity, to the time in- 
terval and a factor (1 - y/k) with k as a value for the 
environment capacity . 

dy = dx • r • y • (1 - y/k) 

r is a prop, factor, the growth rate. 


dy/dx = r y (1 - y/k); r = 0.5; k = 3000; y(0) = 100 
#49: spec_sol_sep 


f 9 1 


1, 0.5 y 1 

„ 0, 100 

. [ 3000 , 



#50: SOLUE |2LN{y) - 2 LN 

x/2 


y - 3000 
29 


+ 2 n - i = x, y 


#51: y = 


3000 e 


x/2 

e +29 


(Initial population) 


I let DERIVE compute the limit of this function to obtain the satu- 
ration quantity 
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dy/dx = r y (1 - y/k) ; r = 0.5; k = 3000; y(0) = 100 (Initial, population) 

V 


#49: spec_sol_sep 1, 0.5 y- i 


#50: S0LUE |2 -LN(y) - 2 LN 

x/2 


3000 

y - 3000 
29 


#51: y = 


#52: lin 


3000 e 


x/2 

e +29 


3000 e 


x/2 


= 3000 


x+m x/2 

e +29 

#53: f e Id ^0 . 5 y |l - ■ 

-30 <= t <= 30 
#54: [[0, 100], 3000] 


3000 


, 0 , 100 


+ 2 - n - i = x, y 


x, y, 0, 20, 2, 100, 2900, 200 



Example 5: The general equation expressing the relation between 

electromotive force U and current I in an circuit con 
taining the resistance R and the inductance L is: 

U = L • I ' + R • I 

Solve this equation first for a constant U = UO. 

I (t = 0) =0 

Then for U 0 = 20V, R = 5 Ohm, L = 0.1 Henry! 

Solve the equation for U = U 0 sin (w t) ! (w = 40) 
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#55: u0 := 

[ u0 - r - i 

1 , , t , i 

1 

1 -LN(i -r - u0> 


f 1 LN(i r - u0> 

#58: SOLUE t = c, i. Real 


r - (c + t ) 

#59: i = IF -n < IM i 

1 


c r/1 - r t/1 


- c r/1 - r t/1 


#60: i = 


[ u0 - r - i 

1, , 0, 0, t, i 

1 

1 LN(-u0) 1 LN(i r - u0) 

#62: = t 


f 1 LN(-u0) 1 LN(i r - u0) 

#63: SOLUE = t, i 


- r t/1 


u0 u0 - e 


#64: i 


20 20 e 


- 5 t/0.1 


#65: i = 


- 50 t 

#66: i = 4 - 4 e 

The next question in this example leads to a linear monic differen- 
tial equation of the form y 1 + p (x) • y = q(x) . The function 

LINEARl (p, q, x, y, xO , yO ) helps us to find the solution. My special 
form LINl (p, q, x, y ) gives the solution containing the parameter c. 


LINEARl(p, q, x, y, x0, y0) := y 


T ^INT(p, x, x0, x) 
y0 + J q-e dx 

X0 

INT(p, x, x0, x) 


r r u0 

#68: LINEARl , SIN(wt), t, i, 0, 0 

, 1 1 


#69: i = 


- r t/1 r t/1 

u0 e (1 w - e (luCOS(tu) - r SIN(t w))) 

2 2 2 
1 -w + r 


- 50 t 

#72: i = 1.95 e - 1.95 COS(40 t) + 2.43 SIN(40 t) 
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[ .I .lx 

J q - e dx + c 

#73: linl (p, q, x, y) : = y = 

J p dx 


f r u0 SIN(u t) 

#74: linl „ „ t, i 

[ 1 1 



r-t/l 

u0 (1 w C0S(t w) - rSIN(t w)) 

#75: 

J. - U C 


2 2 2 
1 -w + r 

#76: 

- 5 0/0.1 

0 = ce 

20 (0.1 40 COS(0 40) - 5 SIN(0 40)) 

2 2 2 
0.1 40 + 5 

#77: 

NSOLUE 

- 5 

0 = ce 

0/0.1 20 (0.1 40 COS (0 40) - 5 SIN{0 40)) 

, c. Real 

2 2 2 
0.1 40 + 5 

#78: 

c = 1.95 


#79: 

- 50 t 

i = 1.95 e 

- 1.95 COS{40 t) + 2.43 SIN{40 t) 


#80: i = - 1.95 COS(40 t) + 2.43 SIN(40 t) 


This is the plot with both solutions #66 (black) and #72 = #79 (red) . 
The exponential term in #79 is called transient, because it usually 
becomes negligibly small after a short lapse of time. Graph of #80 
(blue) shows the remaining steady-state term. 



(Annotating plots is also a benefit of later Derive releases together with direct simplification using the 
= sign at the end of an expression.) 
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Example 6: The ODE y* = x + y defines a nice direction field. 

Let DERIVE plot this field including a family of curves 
fitting into it. 


#81: linl(-l, x) 

x 

#82 : y = ce - x - 1 

#83: DIRECT I 0N_FIELD(x + y, x, -8, 8, 16, y, -6, 6, 12} 

#84: UECT0R(y = c e* - x - 1, c, -10, 10) 

x 

#85: UECT0R(y = c e - x - 1, c, 0, 0.2, 0.02) 
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101/91 Belfast, GBR 

102/91 Montabaur, D 

103/91 Herford, D 

104/91 Stockholm, S 

105/91 Six Fours les Plaques, F 

106/91 Weyhe, D 

107/91 Ludwigshaf en, D 

108/91 Hohenbrunn, D 

109/91 Malm", S 

110/91 Eindhoven, NL 

111/91 Annemasse, F 

112/91 Wien, A 

113/91 Wien, A 

114/91 Grahamstown, SA 


115/91 Milano, I 
116/91 Middlesbrough, GBR 
117/91 Segrate, I 
118/91 Wroclaw, PL 
119/91 Paisley, GBR 
120/91 Graz, A 
121/91 Espoo, SF 
122/91 Waldburg, D 
123/91 London, BR 
124/91 Oberdorf, CH 
125/91 Worms, D 
126/91 Eindhoven, NL 
127/91 Madrid, E 
128/91 Weinheim, D 


If there is any new member, who has not received D-N-L #1, please 
inform the Editor. You will get #1 together with D-N-L #3. 

For the next issue it is planned to present a "List of Members". 
If anybody does not want to be mentioned by name, please inform 
the editor. 



scanned from the original DNL#1 
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VECTOR ( [ -5 • COS ( t) - c-COS(-5t), -5-SIN(t) - c • SIN (-5t) ] , c , -3 , 3) 


